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A systematic numerical study ol non-pairwise vortex interaction forces in the Ginzburg-Landau 
model for single- and multicomponent superconductivity is presented. The interactions are obtained 
by highly accurate numerical free energy minimization. In particular a three-body interaction is 
defined as the difference between the total interaction and sum of pairwise interactions in a system 
of three vortices and such interactions are studied for single and two-component type-1, type-2, and 
type- 1.5 superconductors. In the investigated regimes, the three-body interaction is found to be 
short-range repulsive but long-range attractive in the type-1 case, zero in the critical n (Bogomoln'y) 
case, attractive in the type-2 case and repulsive in the type- 1.5 case. Some systems of four vortices 
are also studied and results indicate that four-body forces are of substantially less significance than 
the three-body interactions. 
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I. INTRODUCTION 

Since the discovery by Onsagcr in supcrfluids and by 
Abrikosov in superconductors^, vortices have been found 
to play crucial role in many aspects of superconductiv- 
ity and superfluidity^. The main difference which can be 
observed in different types of superconductors is how the 
vortices interact and the interaction between a vortex 
pair is known to be attractive for type-1 but repulsive 
for type-2 superconductors^. In a multicomponent su- 
perconducting system it is possible to have also type- 1.5 
behaviour with short-range repulsive but long-range at- 
tractive interactions^^, for a review see^. While type-2 
superconductors exhibit a vortex state with vortices in an 
Abrikosov lattice, type- 1.5 superconductors also exhibit 
a semi-meissner state where the non-monotonic interac- 
tions lead to more complex formations such as vortex 
clusters. 

Non-linearity in the Ginzburg-Landau model implies 
that interactions between many vortices cannot be de- 
scribed simply as the sum of pairwise interactions, al- 
though it is a good approximation for forces between 
widely separated vortices. It can lead to quantitative 
changes in vortex lattices in type-2 superconductors. It 
was sugg ested irP that non-pairwise contributions to in- 
teractions between many vortices in certain type- 1.5 sys- 
tems can dramatically affect vortex configurations. The 
possibility of non-trivial effects caused by non-pairwise 
interactions motivates a further study of such interac- 
tions as done here. In particular the three-body interac- 
tion defined as the difference between total and sum of 
pairwise interactions between three vortices will be stud- 
ied. Also in the more simple single-component type-1 
and type-2 cases, the form of three-body interactions are 
known only for a very limited set of configurations. In 19 
the problem of three vortices in a single-component type- 
2 and type-1 system were studied only for the case of an 
equilateral triangle. In the two-component case, three- 
body intervortex forces were studied for one fixed posi- 



tion of two vortices as a function of arbitrary position of 
a third vortex onljEl Thus multibody intervortex forces 
are known only for few very limited set of configurations 
and it warrants a more extensive study. It will there- 
fore be investigated in Sec. |IV A[ Sec. |IV C| and Sec. 
|IVD| for a wider set of configurations than previously 
studied. We also present the first study of four-body 
forces thereby assessing the scaling of multibody forces 
when extra vortices are included. The known analytical 
suppression of three-body interactions in the Bogomoln'y 
limitPSl, k = will be used as a benchmark to check 
the accuracy of our numerical scheme in Sec. 



systems of four vortices are studied in Sec. |IV 



TVBland 



The free energy density of a two-component supercon- 
ductor in Ginzburg-Landau theory is 



3 = 1,2 

+ 1 (V x A) 2 - t^H^Icos^ - ip x ), (1) 



where q is the electric charge of the superconducting 
charge carriers, ipj = \ipj\e lV: > are the order parameters 
of the two components, A is the magnetic vector po- 
tential and rj is the Josephson coupling, atj and {3j are 
Ginzburg-Landau coefficients corresponding to each com- 
ponent. To determine intervortex forces the free energy 
in Eq. [I] should be minimized with respect to both of the 
order parameters as well as the vector potential. This 
is done here numerically using a finite difference method 
developed irP and briefly described in Sec. [n] Single- 
component systems are described by the same free en- 
ergy except without the terms corresponding to the sec- 
ond condensate and the coupling. Physical realization of 
two-component models lik e the one above are multiband 
superconductors^- 5118 1 21 ' 22 1 and pr ojected superconducting 
states of metallic hydroge n 1 23 * 24 ! 
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II. NUMERICAL METHOD 

The numerical method used to calculate intcrvortex 
forces is essentially that used also irP. A starting guess 
is done for A and ipj, setting the phase winding and 
positions of vortices. The free energy in Eq. [T] is then 
minimized by Newton- Raphson based iteration. Vortex 
pinning is implemented by setting ip = at the position 
of vortices which is minimally invasive, and thus highly 
accurate for determining multibody forces, but does not 
work for too short vortex separations as vortices will then 
tend to escape from their positions while leaving singular 
points with ipj = at the original positions due to the 
pinning constraint. Whether this happens can easily be 
checked by plotting the solutions for A and ipj . 

The energy was minimized for a rectangular grid with 
size Ni until it converged to a value E(Ni). This was 
then repeated for a new grid size N 2 , typically a factor 
four (double in x- and y-direction) greater, to give an 
energy E(N 2 ) etc. This allowed checking for convergence 
with respect to grid size. Final grid sizes were in the 
order N n ~ 10 7 , giving a relative convergence 



E(N n ) - gCjVn-l) 
E{N n ) 



< 10 



-5 



(2) 



Interaction energies are calculated by first finding the 
energy Ei of a system with only a single vortex. The 
pairwise interaction energy E(R) of a vortex pair with 
intervortex distance R is calculated by finding the total 
energy E to t(R) of the system with two vortices and then 
subtracting the single vortex energies so 



E(R) = E tot (R)-2E 1 . 



(3) 



Minimization of the energy leads to the Ginzburg-Landau 
equations 



Similarly, if the total energy in a system of three 
vortices with intervortex distances Ri, R 2 and -R3 is 
Etot(Ri,R2,R3), the total interaction energy of the sys- 
tem is 

Ei n t(Ri, R 2 , ^3) = E to t(Ri, R2, R3) — 3-Ei (4) 

and the three-body interaction is the difference between 
total interaction and sum of pairwise interactions 

E 3 (Ri, R 2 , R3) = E int (Ri, R 2 , R3)—E(Ri)—E[R2)— E(R 3 ). 

(5) 



III. PAIRWISE INTERACTIONS 

For the asymptotic long-range pairwise interaction be- 
tween vortices in a single-band system, analytical re- 
sults are available^!! Introducing ip — \J~^f-.f^ v , Q = 

A — ^ and the Ginzburg-Landau parameter k = y|f^ 
the free energy is 

1 



F= dV 



(1 - f 



K 



Q 2 f + (V x Q) 



K 



I = / (1 - f) 



and 



V x V x Q + / 2 Q = 0. 



(7) 



(8) 



The cylindrically symmetric solutions fo and Qo to Eq. 
[7]and Eq. [8]for a single vortex with boundary constraints 



fo — >• and Q r k as r^O 



(9) 



fo — >• 1 and Q r — > as r -> 00 (10) 
are for large distances r 

Qo = cK^rYe (11) 

and 

g=l-f = dK (V2Kr), (12) 

where Ko and K\ are modified Bessel functions while c 
and d are some values dependent on k. 

In the presence of a second vortex, far away from the 
first one, the solutions to Eq. [7] and Eq. [8] change to 



Q = Q0 + Q1, 



/ = /o + /i, 



(13) 



(14) 



where Qi and /1 are small perturbations. Inserting Eq. 
[13] and Eq. [14] into Eq. [7] and Eq. [8] yields perturba- 
tion equations for Qi and f\ and upon linearization the 
solutions are 



Q 1 (f) = Q (r-r') ! 



/ 1 (r)= 5 (r-r'), 



(15) 



(16) 



where r' is the position of the second vortex. That is, 
the cylindrically symmetric solution for the first vortex 
is perturbed by the asymptotic long-range solution for 
the second vortex. 

To find the interaction energy, Eq. 13 and Eq. 14 are 
inserted into Eq. [6] This yields 



F = U + U int , 
where Uq is the single vortex energy 



Un = dV 



(0) 



(17) 



Qo/o + (V x Qo) 



(18) 
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and f/int is the interaction energy which consists of two 
terms so Ui nt — U1 + U2 where 



Ui = 2 <j>dS ■ [Qi x (V x Q )] 



2 ^(1S- ( /; 



(19) 



(20) 



and integration is done over a surface containing the first 
vortex and far away from both vortices. Evaluation of 
these integrals yields 



Ui(r) = 2TTC 2 K (r) 



U 2 (r) = -—K (V2 K r), 

K 



(21) 



(22) 



where r is the intervortex distance. Here U\ is a repulsive 
electromagnetic interaction and U2 is an attraction due 
to overlapping of the vortex cores. 
It can be showrP that 



c 


for 




1 




K > 




c 


for 




1 




K = 




c 


for 




1 


71 


K < 





(23) 



Hence there is an attraction between type-1 vortices (n < 
^j), repulsion between type-2 vortices (n > and 
no interaction between vortices in the Bogomoln'y limit 



IV. RESULTS 

This section contains results of the numerical calcula- 
tions of three and four-body interactions according to the 
method in Sec. [II] Three-body interactions are presented 
for single-component type-1, critical n and type-2 as well 
as two-component type-1. 5 systems in Sec. |IV A| - [TVD| 
respectively and four-body interactions are presented in 
Sec. |IVE| Three-body interaction energies are presented 
in Fig. [1] Fig. [3] and Fig. [5] as function of the posi- 
tion of the third vortex when a first vortex pair is fixed 
at positions (±^,0). Four-body interactions are shown 
as function of side length R in a square configuration in 
Sec. |IVE| All energies arc given in terms of the energy 
of a single vortex in the system studied. The "minimally 
invasive" vortex pinning mechanism does, as explained in 
Sec. [H] not work for short distances and hence no data is 
available for short vortex separations, which is the reason 
for the empty regions in Fig. [T|- Fig. [6j 



Single-Component Type-1 Three-Body 
Interaction 



Fig. [T] shows the total interaction energy, the sum of 
pairwise interaction energies and the three-body interac- 
tion energy of a single-component type-1 system. The 
third row is hence the difference between the first two 
rows. The parameters used are a = — 1, /3 = 1 and 
q = 2.5. Fig. [2] shows the same data for total interac- 
tion and three-body interaction as that in Fig. [I] but as 
contour plots. 

What can be noted from the results in Fig. [T]is that the 
three-body interaction is repulsive when all vortices are 
close together but also has an attractive region at longer 
vortex separations. The pairwise interaction is however 
attractive at any distance as expected from known ana- 
lytical results. The plots for different values of R\ show 
how the maximum magnitude of the three-body inter- 
action decreases with increasing R\. The magnitude of 
the repulsive part in the three-body interaction appears 
to decrease more strongly with increasing R\ than the 
attractive region does. 



Single-Component Critical Kappa Three-Body 
Interaction 



It can be shown analytically that the pairwise interac- 
tion energy is zero between a vortex pair i n a su percon- 



ductor with k at the critical value n c — JJ 3 | 26 l This is 
consistent with numerical results obtained here by cal- 
culations of interaction energies in a critical k supercon- 
ductor. AlscE^l there is no multibody interaction between 
vortices in a system with critical k. Numerical results 
obtained here suggest that indeed no such three-body in- 
teraction exists in a critical k system, at least within the 
order of 1CP 5 £i where E\ is the energy of a single vortex. 
Thus this gives an estimate of the accuracy of our nu- 
merical scheme. The case studied is a single-component 
system with parameters a = — 1, ft = 1 and q = 2. 



Single-Component Type-2 Three-Body 
Interaction 



Fig. [3] shows the total interaction, sum of pairwise 
interaction and three-body interaction energy of a single- 
component type-2 system with a first vortex pair placed 
at a distance i?i from each other. The interaction energy 
is shown as a function of the position of the third vortex. 
The parameters used here are a = — 1, = 1 and q = 1.5. 
Fig. [4] shows the same data for total interaction and 
three-body interaction as that in Fig. |3]as contour plots. 

In Fig. [3] it is seen that all the type-2 systems stud- 
ied exhibit repulsive pairwise and total interactions con- 
sistent with known analytical and experimental results. 
The three-body interaction on the other hand is at all 
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(g) Three-body interaction. Ri = 1.2. (h) Three-body interaction. i?i = 1.6. (i) Three-body interaction. i?i = 2.0. 

Figure 1: Interactions in a single-component type-1 superconductor with Ginzburg-Landau parameters a = — 1, 
/3 = 1 and q = 2.5. Energy is shown as function of position of the third vortex with a first vortex pair fixed at 
(±^,0) (shown by black lines). Each column corresponds to a certain value of 



distances attractive causing a reduction in the total in- 
teraction energy compared to the sum of pairwise interac- 
tion. This is consistent with three-body forces in type-2 
two-component models^, however, here a much wider set 
of configurations is studied. 



D. Two-Component Type-1. 5 Three-Body 
Interaction 

Multibody forces in the two-component case should in 
general not be expect to be the same as in the single- 



component case. This is because in the two-component 
case a vortex is a composite object made of two frac- 
tional vortices. For a review see^. Thus the multibody 
forces which arise due to mutual deformation of vor- 
tices can include splitting of a composite vortex into con- 
stituent fractional vortices. The fractional vortices have 
a number of exotic properties such as flux delocation and 
inversion^ which makes it an especially complicated non- 
linear problem. Fig. [5] shows the total interaction, sum of 
pairwise interactions and three-body interaction energy 
of a two-component type- 1.5 system with a first vortex 
pair placed at a distance R\ from each other. The inter- 
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action energy is shown as a function of the position of 
the third vortex. The parameters used are ai = —1.0, 
pi = 1.0, a 2 = 3.0, P 2 = 0.5, q = 1.5 and rj = 7.0. Fig. § 
shows the same data for total interaction and three-body 
interaction as that in Fig. [5] but as contour plots. 

All of the type- 1.5 systems studied here show a pair- 
wise and total interaction which is short-range repulsive 
and long-range attractive while the three-body interac- 
tion is at all distances repulsive. These observations are 
consistent with results irP but a wider range of vortex 
configurations have been studied here, giving a complete 
picture of the three-body interaction. In Fig. [5] it can 
be seen how the maximum magnitude of the three-body 
interaction decreases with increasing R\. 



E. Four-Body Interactions 

It should not in general be possible to describe the 
interaction between more than three vortices as neither 
the sum of pairwise interactions nor the sum of pairwisc 
and three-body interactions. It is therefore of interest to 
investigate also the interactions in systems of more vor- 
tices and hence a four-body interaction is defined as the 
difference between the total interaction and the sum of 
pairwise and three-body interactions in a system of four 



vortices. Such interactions are studied here for the same 
single-component type-1 and type-2 and two-component 
type-1.5 systems as studied in the previous sections. As 
a system of four vortices has more degrees of freedom 
than three vortices, the case studied is limited to that of 
a square configuration. 

Fig. [7] shows interaction energies of a configuration of 
four vortices in a square with side length R. The top row 
shows interaction energy for a single pair with distance 
R, a single triangle with sides R, R and \/2R as well 
as the four-body interaction of the square. The bottom 
row shows the total interaction of the square, the sum 
of pairwise interactions as well as the sum of pairwisc 
and three-body interactions. Left column is the same 
single-component type-1 system previously studied with 
parameters a = — 1, /3 = 1 and q — 2.5. Middle column is 
a single-component type-2 system with parameters a = 
— 1, f) = 1 and q = 1.5. Right column shows the two- 
component type-1.5 system with parameters a\ = —1.0, 
fix = 1.0, a 2 = 3.0, /3 2 = 0.5, q = 1.5 and r) = 7.0. 

The most important conclusion to be drawn from Fig. 
[7] is that, in the studied regimes, the four-body inter- 
action is smaller than pairwise and three-body interac- 
tions, except possibly in the type-1 case where in some 
cases it is of similar size as the three-body interaction. 
In the bottom row of Fig. [7] it is seen however that in 



(a) Total interaction. 



(b) Total interaction. 



(c) Total interaction. 




(d) Sum of pairwise interactions. (e) Sum of pairwise interactions. (f) Sum of pairwise interactions. 




y^^x y^-4x y^-4x 

(g) Three-body interaction. Ri = 1.2 (h) Three-body interaction. Ri = 1.6 (i) Three-body interaction. Ri = 2.0 

Figure 3: Same type of data as in Fig. [TJbut for a single component type-2 system with parameters a = — 1, /3 = 1 

and q — 1.5. 



all cases the approximation of the total interaction by 
the sum of pairwise and three-body interactions is signif- 
icantly better compared to that of only the sum of pair- 
wise interactions, especially at short distances. Adding 
also the four-body interaction does not make as big dif- 
ference as the three-body interaction. This would imply 
that the three-body interaction is of greater importance 
than four-body interactions in the systems studied here. 
Even in the type-1 case where the four-body interaction 
seems stronger, it is not of as great importance since in 
a system of four vortices there is only one contribution 
of four-body interaction but four contributions of three- 
body interaction. 



V. CONCLUSIONS 



Non-pairwise interactions between vortices in single 
component type-1 and type-2 as well as two-component 
type-1. 5 superconductors have been studied numerically 
in Ginzburg-Landau theory. In the studied parameter 
sets the results show that there is a three-body interac- 
tion which is short-range repulsive but weakly long-range 
attractive in the type-1 case, attractive in the type-2 case 
and repulsive in the type- 1.5 case. In the critical case, 
k = it is confirmed that just as there is no pair- 
wise interaction, there are also no three-body interac- 
tions (within the accuracy of our numerical scheme) . Re- 
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(d) Three-body interaction. R± = 1.2. (e) Three-body interaction. Ri = 1.6. (f) Three-body interaction. Ri = 2.0. 

Figure 4: Contour plots of the data from Fig. [5] 



suits presented for systems of four vortices imply that the 
three-body interactions are of greater importance than 
four-body interactions so that inclusion of a three-body 
interaction might be an improvement to the approxima- 
tion of the total interaction by the sum of pairwise in- 
teractions. The results are important for giant vortices 
in type-1 superconductors^, pinned clusters of vortices 
in type-2 superconducto rs a nd vortex cluster states in 
type-1. 5 superconductors^ 
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(d) Type-1. (e) Type-2. (f) Type-1.5. 

Figure 7: Interaction energies for four vortices in a square with side R. Total interaction, sum of pairwise 
interactions and sum of pairwise and three-body interactions, (a)-(b) show single-component type-1 and type-2 
systems with parameters a — —1, /3 = 1 and q — 2.5 and a = —1, f3 — 1 and q — 1.5 respectively, (c) is a type-1.5 
system with parameters ot\ — —1.0, /?i = 1.0, a 2 — 3.0, (3 2 — 0.5, q — 1.5 and r\ = 7.0. 



